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The origin of an extended X-ray emission apparently associated with the globular 

cluster 47 Tucanae 
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Abstract 

Using the Suzaku X-ray Imaging Spectrometer, we performed a 130 ks observation of an extended X-ray emis- 
sion, which was shown by ROSAT and Chandra observations to apparently associate with the globular cluster 47 
Tucanae. The obtained 0.5 — 6 keV spectrum was successfully fitted with a redshifted thin thermal plasma emission 
model whose temperature and redshift are 2.21q 3 keV (at the rest frame) and 0.34 ± 0.02, respectively. Derived 
parameters, including the temperature, redshift, and luminosity, indicate that the extended X-ray source is a back- 
ground cluster of galaxies, and its projected location falls, by chance, on the direction of the proper motion of 47 
Tucanae. 
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1. Introduction 



hard that they require power-law models with photon indices 
of r = 2.1 ± 0.3 (47 Tuc) and T 2.0 ± 0.2 (NGC 6752), 
or a thermal plasma emission model with temperatures of 
Galactic halo with a typical velocity of - 200 km s^^, fcT = 3.7ti:3 keV (47 Tuc) and fcT^ 2.9l:o:? keV (NGC 6752) 



Many globular clusters in the Galaxy move through the 



which exceeds the sound velocity (a few tens of km s to 
150 km s^^) specified by roughly estimated halo plasma tem- 
peratures (T - 0.7 - 1.4 X 10^ K by Savage & de Boer 1981; 
T - 1.5 - 1.6 X 10^ K by Pietz et al. 1998). Then, a bow shock 
is expected to form between the halo plasma and any gas (intra- 
cluster gas) in a moving globular cluster (Ruderman & Spiegel 
1971). Since the temperature of the post-shock plasma should 
then become ~ 10^ K (e.g., Okada et al. 2007), we expect 
to detect diffuse X-ray emission with a shape that traces the 
bow shock. Several X-ray satellites have been observing glob- 
ular clusters in search for such bow-shock-heated X-ray emit- 
ting plasmas. Indeed, using the Einstein satellite, Hartwick, 
Cowley, & Grindlay (1982) first detected such diffuse emis- 
sions around 47 Tucanae (hereafter 47 Tuc), lo Centauri, and 
M22, although a part of the emissions was resolved into point 
sources by later observations (Koch-Miramond & Auriere 
1987; Krockenberger & Grindlay 1995; Gendre, Barret, & 
Webb 2003). Subsequently, Krockenberger & Grindlay 
(1995) newly reported diffuse emissions in 47 Tuc, followed 
by, possible detections of such a diffuse emission in several 
globular clusters; e.g., Hopwood et al. (2000) in NGC 6779, 
and Okada et al. (2007) in 47 Tuc, NGC 6752, and M5. 

Among the extended emissions so far detected, those in the 
globular clusters 47 Tuc and NGC 6752 are of particular inter- 
est. They spatially coincide with the directions of projected 
proper motions of these globular clusters (Krockenberger & 
Grindlay 1995; Okada et al. 2007), and were hence consid- 
ered to have physical relationships with the globular clusters. 
According to Okada (2005) and Okada et al. (2007), the 0.5- 
4.5 keV Chandra spectra from these extended sources are so 
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which largely exceed values expected from the bow shock heat- 
ing (^ 10^ K). While the two X-ray sources have no optical 
identifications (Krockenberger & Grindlay 1995; Okada et al. 
2007), Okada et al. (2007) reported that both have possible ra- 
dio counterparts in the 843 MHz Sydney University Molongo 
Sky Survey (SUMSS; Bock et al. 1999). 

From these properties, the extended emissions apparently as- 
sociated with 47 Tuc and NGC 6752 were considered to arise 
via inverse Compton scattering (Krockenberger & Grindlay 
1995) or non-thermal bremsstrahlung (Okada et al. 2007) of 
high energy {E 2Q— 100 keV) electrons that are stochasti- 
cally accelerated in the bow shock. The interpretation is attrac- 
tive because the shock is expected to be a moderate one with a 
Mach value of 10, and the condition is much different from 
those in the more typical acceleration sites such as supernova 
remnants and jets of active galactic nuclei. 

As an alternative explanation to those extended X-ray 
sources, Krockenberger & Grindlay (1995) and Okada et 
al. (2007) also considered a chance coincidence with a back- 
ground cluster of galaxies that is not related to the globular 
clusters. This alternative must be kept in mind, even though 
its possibility was estimated low (< 1% by Krockenberger & 
Grindlay 1995 and Okada et al. 2007) in 47 Tuc. 

In the previous spectral analysis of the Chandra data from 
the extended emission in 47 Tuc and NGC 6752, Okada et al. 
(2007) were unable to distinguish a power-law model from a 
thermal emission model because of rather large statistical er- 
rors. In the present paper, we utilize the larger effective area 
and lower background level of Suzaku, to perform detailed 
spectral analysis of the extended emission of 47 Tuc. Based 
on the model fitting result, we conclude that the emission is 
from a background galaxy cluster with a redshift of 0.3, and 
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the spatial coincidence between the extended emission and the 
globular cluster is accidental. 

Throughout the paper, cosmological parameters of flu = 
0.28 and Hq ~ 70 km s~^ Mpc"^ are used in calculations. 

2. Observation and Data Reduction 

The globular cluster 47 Tucanae was observed with Suzaku 
(Mitsuda et al. 2007) on 2007 June 10-12 (observation ID 
502048010). Since the target of this observation is the extended 
emission (EE) at the north east region of 47 Tuc, we set the 
nominal pointing of the satellite at (00'^24'"50^ -71°59'46"), 
which is ^ 6' off the cluster center. 

In the present study, we focus on the data taken with the X- 
ray Imaging Spectrometer (XIS; Koyama et al. 2007), which 
comprises four X-ray charge coupled device (CCD) sensors 
each placed on the focal plane of the X-ray Telescope (XRT; 
Serlemitsos et al. 2007). The four pairs of XIS and XRT are 
co-aligned together, and have the same field of view (FOV) 
of 17.'8 X 17'.8. Since one of three front-side illuminated (PI) 
CCD chips, XIS2, was not operational since 2006 November, 
we utilized the data from the remaining two FI sensors (XISO 
and 3) and a back-side illuminated (BI) one (XISl). In the 
present observation, the XIS was operated in the normal mode 
without any window or burst option, but incorporating the 
spaced-row charge injection method (Nakajima et al. 2008) to 
restore the energy resolution of the CCDs. 

After removing periods of the Earth elevation angle less 
than 5° (ELV< 5°), the day Earth elevation angle less than 
20° (DYE_ELV< 20°), and the South Atlantic Anomaly, we 
achieved a net exposure of 132 ks. Flickering pixels were re- 
moved from the data by using cleansis version 1.7. Then, 
cleaned event files were generated employing the same event 
extraction criteria as in the Suzaku pipe line processing (ver- 
sion 2). 

The present data reduction and analysis were performed 
using HEADAS package version 6.4.1 and XSPEC version 
11.3.2. In spectral fitting, redistribution matrix files and an- 
cillary response files (ARFs) for the XIS/XRT were generated 
using xisrmfgen version 2007-05-14 and xissimarf gen 
(Ishisaki et al. 2007) version 2008-04-05, respectively, with the 
calibration files which are provided by the calibration database 
(CALDB) version 2008-04-01. In the spectral fitting described 
below, we ignored data in the 1.8 — 1.9 keV band so as to avoid 
calibration uncertainties around the Si-K edge. 

Events with energies above 10 keV, taken with the Hard X- 
ray Detector (HXD; Takahashi et al. 2007), were not utilized in 
the present analysis. This is because the HXD lacks imaging 
capability, and hence we cannot exclude contamination by X- 
rays from a number of point sources associated with 47 Tuc 
(eg. Verbunt & Hasinger 1998; Grindlay et al. 2001; Heinke et 
al. 2005). 

Since the spatial resolution of the XIS/XRT is ~ 2', we can- 
not exclude, using the XIS data alone. X-ray point sources that 
overlap with the EE. To determine their spectral shapes and 
fluxes, we also utilized archived Chandra ACIS data of 47 Tuc 
acquired in 2000 March for a total exposure of about 70 ks 
(obsid=953 and 955). We used CI AO (Chandra Interactive 
Analysis of Observations) version 4.0.2 and CALDB version 



3.4.5 to extract point source spectra. Like in the Suzaku data 
analysis, we also used XSPEC when performing model fitting 
to the ACIS spectra. 

3. Image Analysis 

3.1. Soft and hard band images 

In Fgure 1, we present images obtained with the FI cam- 
eras (XISO and 3) in the soft (0.5 - 1.5 keV) and hard (1.5 - 
6.0 keV) bands, after subtracting the non X-ray background 
(NXB) and correcting for vignetting and exposure. We esti- 
mated the NXB of the XIS using dark (night) Earth data taken 
within ±150 days of our observation of 47 Tuc. The night 
Earth data were summed up, with weights according to geo- 
magnetic cut-off rigidity which the spacecraft experienced at 
the data acquisition. This was performed by xisnxbgen 
(Tawa et al. 2008). Then, we created NXB images in the 
soft and hard bands, and subtracted them from the raw im- 
ages. After subtracting the NXB, we smoothed each image 
with a two-dimensional Gaussian of ct = 6". The diffuse X- 
ray backgrounds, namely the cosmic X-ray background (CXB) 
and Galactic diffuse emission, are still included in the images. 

In figure 1, we see several point sources, and the very bright 
47 Tuc core region which consists of numerous X-ray point 
sources (e.g., Heinke et al. 2005). At the center of the two 
images, we also observe a clear concentration of X-ray events 
as Krockenberger & Grindlay (1995) and Okada et al. (2007) 
reported. Thus, we reconfirm the EE phenomenon with the 
Suzaku data. To extract photons from the EE region, we define 
a circular region (white circle in figure 1) with a radius of 150", 
centered on (00'^24"i44?2, -71°59'33."5) where Okada et al. 
(2007) found the maximal surface brightness. 

As indicated with a black solid circle and a label "PS" in 
figure 1, a faint point source is recognized at the north west rim 
of the event extracting region. Although the EE is still apparent 
in the hard band image, the point source is no longer visible 
therein. At a consistent position, we find a point source also in 
the ACIS image. Therefore, we consider that the XIS source 
is not a brightness enhancement associated with the EE, and 
hereafter exclude it using a circular region with a radius of 1'. 
This region is expected to enclose 50% of X-ray photons from 
the point source, while the remaining half will fall out of the 
region; a half of those photons (25% of the total flux from the 
source) are in turn estimated to fall inside the event extracting 
region around the EE, and contaminate the EE spectrum. This 
effect is considered later in section 4. 

3.2. Radial profile of the extended emission 

In the previous studies by Krockenberger & Grindlay (1995) 
and Okada et al. (2007), the EE was concluded to be extended 
over an angular radius of ~ 2' (2.7 pc assuming a 4.6 kpc dis- 
tance to 47 Tuc). To examine the spatial extent of the emission, 
we calculated its azimuthly averaged radial profile using the 
NXB-subtracted and vignetting-corrected 0.5 — 6.0 keV XIS 
FI image. This was done utilizing a series of annular extracting 
regions, each with 30" width, which are concentric with the 
original event extraction region. The result is shown in figure 2 
after subtracting the CXB and Galactic diffuse background rate 
of 9.4 X 10^^ counts s^^ pixeP^, which we estimated using 
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Fig. 1. Soft (0.5 - 1.5 keV; left) and hard (1.5 - 6.0 keV; right) band images of 47 Tuc taken with XIS FI (XISO plus XIS3), shown after removing the 
two corner regions irradiated with the calibration source. The images are scaled in units of 4 X 10~^counts s"'^ pixel"^. The non X-ray background 
is subtracted using the night Earth image (see text), followed by vignetting and exposure correction, although the diffuse X-ray background is included. 
The white circle with a radius of 150" is the event extracting region for the extended emission (EE), and the black small one shows a soft point source 
(PS) which is excluded from the spectral analysis. 
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Fig. 2. Azimuthly averaged radial profiles of the EE (filled circles) 
shown in units of counts s^^ pixel^^, after subtracting the CXB and 
Galactic diffuse components (see text). The PSF of XIS FI (open 
squares) is also plotted, being normalized to have the same maximum 
value as that of the EE at the innermost annulus. 

another region of the CCDs with no evident point sources. 

In the same figure, we also plot a radial profile of the point 
spread function (PSF) of the XIS/XRT, calculated at the center 
of the EE, and averaged over XISO and XIS3. Thus, the EE is 
clearly more extended than the PSF even though the latter is 
much broader than those of ROS AT and Chandra. 

4. Spectral Analysis 

4. 1 . Extraction of spectra 

We extracted XIS FI and BI spectra of the EE using the 
event extracting region shown in figure 1 (the white circular 
region but excluding the black circle). The 0.5 — 6 keV band 
count rates measured with XIS FI and BI are 28.3 ± 0.05 x 



10^^ counts s^^ and 22.7 ± 0.04 x 10^^ counts s^^, respec- 
tively, with la statistical errors. 

In addition to the EE which is the subject of the present 
analysis, the spectra also contain events from the NXB, the 
Galactic and extragalactic X-ray backgrounds (altogether, dif- 
fuse X-ray background), and X-ray events from several con- 
taminating X-ray point sources that cannot be resolved with 
the XIS/XRT spatial resolution. We assume the Galactic dif- 
fuse emission to have a uniform brightness across the XIS FOV. 
Although its hard component (Worrall et al. 1982; Koyama et 
al. 1986; Ebisawa et al. 2001; Revnivtsev et al. 2006; Krivonos 
et al. 2007) has a strong concentration toward the Galactic 
plane (with a scale height of l.°5 — 3° ; e.g., Revnivtsev et al. 
2006, Krivonos et al. 2007), and hence a steep brightness gra- 
dient, it can be neglected at this high Galactic latitude of ^ 45° 
of 47 Tuc. 

4.2. Background Spectra 

In the following subsections, we estimate the spectral shapes 
and fluxes of these individual background components, and 
create XIS spectral data for each component. By summing all 
these components, we construct background spectra, and then 
subtract them from the raw XIS spectra of the EE. 
4.2.1. The non X-ray Background 

We derived the NXB spectrum from the same stacked night- 
Earth data as described in section 3.1. Since this component 
depends on the CCD location, we examined spectral differ- 
ences among several circular extracting regions on the night- 
Earth image, with the radius ranging from 150" to 300". Each 
region is concentric with the EE extracting region (white cir- 
cle in figure 1). The derived 0.5 — 10 keV NXB spectra were 
consistent with one another within ~ 3%. Therefore, to min- 
imize the statistical errors of the estimated NXB, we adopted 
the largest extracting region (300" radius) for both XIS FI and 
BI. 
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The constructed NXB spectrum is shown in figure 6 in green. 
The count rate has been scaled to the ratio (~ 4.8) of the NXB 
and signal extracting areas. The 0.5 — 6 keV band count rates 
with la statistical errors are 4.4 ± 0.1 x 10"'^ counts s^^ (XIS 
FI) and 5.3 ± 0.1 x lO^^ counts s'^ (XIS BI). 
4.2.2. The diffuse X-ray background 

The diffuse X-ray background consists of two components; 
the Galactic and the extragalactic emissions. The former com- 
ponent is thought to originate from the Galactic halo and the 
Local Hot Bubble (eg. Cox & Reynolds 1987), and expected 
to appear at energies below ~ 2 keV with its surface bright- 
ness depending considerably on the sky direction. The latter, 
the extragalactic component, has been understood as a super- 
position of numerous extragalactic active Galactic nuclei. The 
spectrum is known to be expressed by a power-law model with 
a photon index of F = 1.4 (e.g. Parmar et al. 1999; Lumb et al. 
2002; Kushino et al. 2002) at least over the 2-10 keV band. 

In order to determine the local diffuse X-ray background in 
the present XIS FOV, we extracted another set of XIS FI and 
BI spectra from the same observation data set of 47 Tuc, but 
applying a mask which excludes point sources, the core region 
of 47 Tuc, and the EE itself. The masked image of XIS FI 
is shown in figure 3, and the NXB-subtracted (as described in 
section 4.2. 1) spectra of the diffuse X-ray background are plot- 
ted in figure 4. We fitted these spectra jointly with a model 
which consists of three diffuse X-ray background components; 
a thermal emission from the Local Hot Bubble plasma (mekal 
model in XSPEC; Liedahl et al. 1995); a thermal emission from 
the Galactic halo plasma (mekal); and a power-law model 
with a fixed photon index of 1.4 to account for the extragalac- 
tic component (powerlaw). The latter two components were 
assumed to suffer the line-of-sight Galactic absorption, with 
the absorbing column density fixed at 5 x 10^" atoms cm~'^ 
(Dickey & Lockman 1990) which is a typical value toward 
the present field. The photoelectric absorption coefficient by 
Morrison & McCammon (1983), wabs model in XSPEC, 
was employed. We assumed that the diffuse background has 
a uniform surface brightness over the XIS FOV, and utilized an 
ARE which was calculated using xissimarfgen with the 
UNIFORM option and r_max= 20'. We left free the temper- 
atures, metal abundances, and surface brightnesses of the two 
thermal models, as well as the photon index and surface bright- 
ness of the power-law model. The model gave an acceptable fit 
with xi = 1-10 {i^ = 175); the best fit parameters are listed in 
table 1. 

As to the power-law component, the best fit model gave 
the 2-10 keV flux of 4.4 x 10"* erg cm~^ s"^ sr"\ 
which is ~ 20% lower than the previously reported values; 
5.4 X 10"® erg em"^ s"^ sr"^ by Lumb et al. (2002), and 
5.7 X 10"® erg cm"^ s"^ sr~i by Kushino et al. (2002). The 
deviation can reasonably be explained by the spatial fluctua- 
tion of the extragalactic emission which can vary by about 20% 
(Kushino et al. 2002) depending on XIS pointings. 

We might directly subtract the diffuse background spectrum 
of figure 4 from that of the EE region. However, this introduces 
some systematic bias because the energy-dependent vignetting 
effect of the XRT (figure 1 1 of Serlemitsos et al. 2007) wiU 
cause not only the observed background brightness but also its 
spectral shape to differ between the two regions; in the present 
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Fig. 3. The XIS FI image after filtering out point sources, the 47 Tuc 
core region, and the EE. The image is not corrected for the vignetting 
or exposure. The events plotted in the image were used in the modeling 
of the diffuse X-ray background in the field of 47 Tuc. 
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Fig. 4. The NXB-subtracted diffuse X-ray background spectra of XIS 
FI (black) and BI (red), extracted from the image in figure 3. The solid 
lines represent the best fit model, while their components are individu- 
ally plotted in dashed (thermal), dot-dashed (thermal with absoiption), 
and dotted (power law with absorption) lines. 

case, the two ARFs, one for the masked region (figure 3) while 
the other for the EE, differ by 10 — 20% (due to energy depen- 
dent vignetting effect) in the 2 — 6 keV range if we normalize 
them at 2 keV. Hence, to avoid this problem, we simulated the 
expected contribution of the diffuse background to the EE ex- 
tracting region using the best fit model explained above and the 
corresponding ARE. In producing the fake spectra, we assumed 
a sufficiently long exposure (10^ s), to suppress statistical er- 
rors. This is allowed because the background components are 
understood from previous observations. The 0.5 — 6 keV band 
count rates of the faked spectra are 8.3 x 10"'^ counts s"^ (XIS 
FI) and 5.9 x 10"^ counts s"i (XIS BI). 
4.2.3. Contamination from point sources 

In the previous study using Chandra (Okada et al. 2007), 
six faint X-ray point sources were found within 2.'5 of the 
EE region. In table 2, we list their positions. Although they 
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Table 1. The best fit model parameters for the diffuse background spec- 
tra. 



Model Parameter Value 



Thermal 1 


kT 


„ ^ „ XTTTTT — T — TT — 

0.17 ±^^:^^ keV 




Z* 


i±+°-9 

^■^=■=-0.03 




St 


i.U(in_Q gg 


Absorption 


iVn* 


5 (fixed) 


Thermal 2 


kT 


0.78+° i^^ keV 




Z* 


03+0 04 

^■"«J-0.03 




Et 


n 7C+2.28 
' °-0.35 


Power law 


r 


1.4 (fixed) 




St 


n OO+0.03 



^ "■'^"-0.03 

Xt 1. 10(175) 



'Abundance in terms of the solar value (Anders & Grevesse 
1989). 

txhe 0.5 — 6 keV band model surface brightness in units of 
10~* erg cm~^ sr~^. Absorption is not connected. 
■fLine-of-sight hydrogen column density in units of 
10^0 cm-2. 
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Fig. 5. The XIS FI spectrum (black crosses) and the best fit power-law 
model (solid line) of the soft point source at the north west of the EE. 
The NXB and diffuse X-ray background ai'e subtracted. Data from XIS 
BI are excluded from the plot for clarity, although they were incoipo- 
rated in the fitting. 



Table 2. The coordinates of the contaminating point sources. 



Source # Coordinate 



1 


(00'^24'°14?51, 


-71°58'50"4) 


2 


(00^24'°38?71, 


-72°00'46"3) 


3 


(00^^24'"34?83, 


-72°00'40"2) 


4 


(00^24'°30?26, 


-72°00'33"8) 


5 


(00*^25'°00?70, 


-71°59'59"9) 


6 


(00^24'°42f63, 


-71°59'22"3) 



were successfully removed in the Chandra case, we cannot do 
so from the present XIS data except for the brightest one de- 
scribed in section 3, because of the broader PSF of the XRT 
than that of Chandra. Therefore, we must model and subtract 
their contributions, like the diffuse background. As explained 
below, we estimate the contribution from the soft point source 
(Source 1; figure 1) using the Suzaku XIS data themselves, and 
those of the remaining five point sources (Source 2-5) using the 
Chandra ACIS data assuming that they are not variable. 

Figure 5 shows XIS FI spectrum of Source 1, extracted from 
the black circle (figure 1), shown after subtracting the NXB and 
the diffuse X-ray background. The FI and BI spectra were fitted 
with an absorbed single power-law model in the 0.5 — 5 keV 
band. As listed in table 3, this gave an acceptable fit with = 
1.10 {v ~ 30). The summed spectrum of the remaining 5 point 
sources was extracted from the ACIS data (section 2), using 
circular regions each 5" in radius. The NXB was extracted 
from another region of the same ACIS CCD with no evident 
point sources. Then, we fitted the summed spectrum with a 
single power-law model in the 0.8 — 6 keV band. The best-fit 
ix^ = 2.54 and u — A) model gave a null hypothesis probability 
of 0.041, and the parameters as listed in table 3. 

To obtain a summed contribution of all the point sources to 
the EE, we then faked the summed spectrum of the 5 point 
sources and the soft source separately, by applying appropri- 



Table 3. The best fit model parameters for the contaminating soft point 
source and five faint sources. 



Model 


Parameter 


Source 1 


Source 2 — 5 


Absorption 
Power law 


r 

fluxt 


9 9+2.3 
^•^-1.3 


(fixed) 
1.7±0.2 
3.2 ±0.3 




xi 


1.10(30) 


2.54 (4) 



*Hydrogen column density in units of 10^^ cm . 
tThe 0.5 — 6 keV band model flux in units of 
10~^* erg cm~^ s~^. Not corrected for the absorp- 
tion. 



ate ARFs to the best fit models described above. In calculating 
the ARF for the 5 point sources, we took an average of indi- 
vidual ARFs weighted by their 0.5 — 5 keV ACIS count rates. 
The ARF for Source 1 was calculated referring to the XIS/XRT 
effective area for X-ray photons which leak into the EE event 
extracting region; the source position was set to be that of the 
soft source (the first row in table 2), whilst it is located out- 
side the EE region (white circle in figure 1)'. Based on the 
faked spectrum, the implied 0.5 — 6 keV band count rates are 
2.9 X 10-3 counts s"! (XIS FI) and 2.3 x 10"^ counts s^^ 
(XIS BI). 

4.3. Model fitting to the Extended Emission Spectrum 

Figure 6 shows the raw EE spectra, in comparison with the 
background components estimated so far Table 4 summerizes 
the estimated 0.5 — 6 keV count rate of each background com- 
ponent. As a whole, the background amounts to about 50% of 
the raw counts in each detector In figure 6, cyan data points in- 
dicate the EE spectra obtained after subtracting the three back- 
ground components. Below, we fit them with several models 

' A ratio of the number of photons which leak into the EE region to that 
of photons falling inside the Source 1 region (black circle in figure 1) is 
23% in the 0.5-5 keV band, which is close to with the rough estimation 
(~ 25%) in section 3.1. 
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Fig. 6. The raw (black) and the background-subtracted (cyan) spec- 
tra of the EE obtained with XIS FI (panel a) and BI (panel b). The 
long-accumulated non X-ray background and the faked diffuse X-ray 
background are plotted in green and red respectively. The blue line 
represents the simulated contamination from the six point sources. 



Table 4. The 0.5 — 6 keV count rates of individual spectral compo- 
nents. 



Component 


Rate (10^^ 


counts s 




XISFI 


XIS BI 


Raw 


28.3 


22.7 


NXB 


4.4 


5.3 


DXB* 


8.3 


5.9 


PS^ 


2.9 


2.3 


BGD* 


15.6 


13.5 


EE§ 


12.7 (45%) 


9.2 (40%) 



'Diffuse X-ray background, 
tsix point sources. 

t Sum of the NXB, diffuse X-ray background, and six 
point sources. 

§ Derived from Raw— BGD. Ratios to the Raw count 
rates are also shown. 
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Fig. 7. Spectral fitting to the XIS FI (black) and BI (red) spectra of the 
EE, with (a) a power-law and (b) a single-temperature thermal models. 

which give different physical interpretations. An ARF for the 
EE was calculated assuming a uniform circular emitting region 
with a radius of 50" based on the Chandra ACIS imaging re- 
sult (Okada 2005; Okada et al. 2007). The FI and BI spectra 
were fitted simultaneously, with the overall model normaliza- 
tion fixed to be the same between them. 

First, we fitted the spectra with a single power-law model 
and a single temperature optically-thin thermal model (apec 
in XSPEC; Smith et al. 2001), each subjected to the interstel- 
lar absorption (wabs) as Okada et al. (2007) did. The fitting 
results are shown in figure 7 and listed in table 5. However, 
neither the power-law nor optically-thin thermal model repro- 
duced the spectra well, with xt = 1-31 {i' = 100) and 1.33 (i^ = 
99) respectively. The obtained photon index (F = 2.9 ± 0.2) or 
the plasma temperature {kT = 1.7 ± 0.3 keV) implies a con- 
siderably softer spectral shape than the previous report (Okada 
2005; r = 2.1 ± 0.3 or kT = S.7±f l keV). In section 5, we 
discuss this discrepancy. 

In figure 7, we notice some spectral structures around 
0.85 keV, 1.5 keV, and 5 keV that cannot be explained by 
the employed models. Suspecting that these structures origi- 
nate from redshifted atomic emission lines, we next fitted the 
spectra with a thermal model that has a free redshift z. The 
fit was improved significantly to xj^ = 1.10 (i^ = 98), and 
yielded the metal abundance and redshift of 0.381q^3 times 
solar and z = 0.34 ± 0.02, respectively. Especially the spec- 
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Energy (keV) 

Fig. 8. The same EE spectra as presented in figure 7, fitted witii a red- 
sliifted thermal emission model. 

tral features at ~ 0.85 keV, ^ 1.5 keV and ~ 5 keV have 
been reproduced successfully by redshifted Fe-L, Si-K and 
Fe-K lines, respectively. Incorporating z thus determined, 
the observed flux can be converted to the intrinsic luminos- 
ity of Lo.5-6keV = 5.5 X 10''^ erg s"\ L2_io keV = 2.8 x 
lO'^^ erg s~^, and Lo.1-200 koV = 1-0 x 10"'^ erg in the 
0.5-6 keV, 2-10 keV, and 0. 1 - 200 keV band respectively. 

Since the XIS background spectrum contains Ka emission 
line from aluminum used in, for example, the XIS housing and 
substrate of the CCD, the line feature at 1.5 keV could be due 
to residual Al-K lines caused by a wrong NXB subtraction. 
To examine this possibility, we also tried spectral fittings with 
NXB spectra rescaled by 5 — 10%. However, the feature at 
~ 1 .5 keV can be seen even after subtracting an NXB spectrum 
that is rescaled up by +10%. Since the systematic error (or 
reproducibility) of the NXB estimation is reported to be 5% 
(Tawa et al. 2008), we consider that the structure to be real 
rather than instrumental. For reference, the fit results remain 
unchanged within the errors even if we ignore the 1.4—1.6 ke V 
range in the fitting. 

In figure 8, we notice fitting residuals both in the XIS FI and 
BI spectra at 3.5 keV. However, they have no corresponding 
background features (figure 6) or redshifted major atomic lines. 
No such features are present in the Chandra spectrum, either 
(Okada et al. 2007). They are hence considered as statistical 
fluctuations. 

The above results strongly suggest that the EE is an extra- 
galactic object, rather than a source associated with 47 Tuc. 
Further considering the extended nature and the thermal spec- 
trum, it is most likely a background cluster of galaxies at 
z — 0.34. In the following section, we examine the galaxy 
cluster interpretation of the EE based on the Suzaku results. 

5. Discussion 

5.7. kT — Lx relation and the counterpart in other wave- 
length 

As shown so far, the spectra of the EE are well described 
by thermal plasma emission with a rest-frame temperature of 
kT — 2.2 keV and a redshift of z = 0.34. Furthermore, as plot- 
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Fig. 9. A kT — Lx relation of clusters of galaxies. Crosses represent 
temperatures and bolometric luminosities of individual galaxy clus- 
ters determined by ASCA observations (data taken from Fukazawa et 
al. 2004). Luminosities were obtained by integrating fluxes over the 
0.1 - 200 keV band. Filled triangle shows the EE of 47 Tuc. 



ted in figure 9, its luminosity and temperature are consistent 
with the known temperature-luminosity relation {kT—Lx rela- 
tion) of cluster of galaxies. Therefore, the EE is most naturally 
interpreted as a background cluster of galaxies at a moderate 
redshift. 

We find no counterpart in the optical (Digital Sky Survey; 
e.g. McLean et al. 2000) or near infrared (Two Micron All 
Sky Survey; Skrutskie et al. 2006) surveys. Using Chandra 
deep survey data, Boschin (2002) however reported more than 
20 candidates of clusters of galaxies that have no optical coun- 
terpart. The present background galaxy cluster is perhaps a 
member of those clusters. A deeper optical imagery will reveal 
the expected galaxy clustering. 

5.2. Comparison with previous reports on the EE 

The galaxy cluster interpretation has been examined as an 
origin of the EE in previous papers as well. Using the logiV — 
log S relation of galaxy clusters, Krockenberger & Grindlay 
(1995) and Okada et al. (2007) estimated probabilities of a 
chance coincidence of the EE and a background galaxy cluster 
emission to be less than 0.5% and 0.6%, respectively. Based on 
such low probabilities, these authors argued that the EE cannot 
be a background galaxy cluster. 

In addition to the above probability estimation, Okada et al. 
(2007) used the following argument to rule out the background 
cluster interpretation. First, they determined the EE tempera- 
ture as kT = 3.7 keV from the Chandra ACIS spectrum. They 
hence assigned a luminosity of Lx = 1.1 x 10^"* erg to this 
putative cluster, using the kT — Lx relation of clusters (e.g. 
Ikebe et al. 2002; Fukazawa et al. 2004). Comparing this Lx 
with the measured flux, the source redshift was estimated as 
z > 0.5, and hence the observed angular core radius of the EE, 

~ 0.'6, was converted to a physical size of > 360 kpc. 
Finally, they concluded this Tc to be too large for a cluster 

In the present study, the use of the Suzaku XIS has en- 
abled us to achieve two major improvements (or revisions) over 
Okada et al. (2007). One is that we clearly detected redshifted 
emission lines, which indicate z ~ 0.34 ± 0.02; the Chandra 
data gave no constraint on z. The other is that we measured 
a significantly lower temperature, kT =1.7 keV if assuming 
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Table 5. The best fit parameters of the EE spectra. 



Model 




r 






Ell 


xl M 


Power law 
Theraml 
Redshifted thermal 


2OI5 

q+4.4 
'^■'-'-3.3 

o.y_4 8 


2.9±0.2 


1.7±0.3 

9 9+0.2 
^■^-0.3 




'•O_0.9 
'•3-1.3 

0.34 ±0.02 7.6tli 


1.31 (100) 
1.33 (99) 
1.10(98) 



*Line-of-sight hydrogen column density in units of 10^'' cm^^. 
t Thermal plasma temperature in units of keV. 
t Abundance in terms of the solar value. 
§Redshift. 

'The 0.5 — 6 keV band model surface brightness in units of 10~^ erg cm~^ sr 



z = 0, or kT ~ 2.2 keV at the rest frame if adopting z = 0.34; 
the latter now satisfies the kT — Lx relation of clusters of 
galaxies (figure 9). In addition, using the redshift, the phys- 
ical core radius is now calculated to be ~ 160 kpc, which is 
reasonable for galaxy clusters. 

As reviewed so far, the difference of our conclusion from 
that of Okada et al. (2007) comes mainly from the discrepant 
EE temperatures, kT — 1.7 ± 0.3 keV measured with Suzaku 
(without correction for the redshift) and kT = 3.7^^ 3 keV with 
the Chandra ACIS. Possible causes of this difference include 
an over estimation of the temperature with Chandra, or an un- 
der estimation with Suzaku, or both. As the former possibility, 
the most likely cause is systematic errors in the NXB subtrac- 
tion. As the latter possibility, we may presume that during the 
Suzaku observation, some soft sources became brighter than in 
the Chandra observation. 

Although the Suzaku data could thus be under-estimating 
the EE temperature, significantly higher values of kT would 
be still consistent with the kT — Lx relation. Furthermore, the 
value of z = 0.34 is not affected, since it is determined by the 
redshifted atomic emission lines. We conclude that the close 
spatial coincidence between the EE and 47 Tuc is accidental, 
and they are not physically associated with each other. 
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